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^; Abstract 

It has been shown that the orbits of motion for a wide class of non-relativistic Hamil- 
tonian systems can be described as geodesic flows on a manifold and an associated dual. 
This method can be applied to a four dimensional manifold of orbits in spacetime associ- 
ated with a relativistic system. We show that a relativistic Hamiltonian which generates 
Einstein geodesies, with the addition of a world scalar field, can be put into correspon- 
dence with another Hamiltonian with conformally modified metric. Such a construction 
could account for part of the requirements of Bekenstein for achieving the MOND theory 
of Milgrom in the post- Newtonian limit. The constraints on the MOND theory imposed 
by the galactic rotation curves, through this correspondence, would then imply constraints 
on the structure of the world scalar field. We then use the fact that a Hamiltonian with 
vector gauge fields results, through such a conformal map, in a Kaluza-Klein type theory, 
and indicate how the TeVeS structure can be put into this framework. 
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1. Introduction 

The Hamiltonian [11 



P. : k = —g^p\ (i) 

> . 2m 

with Hamilton equations (written in terms of derivatives with respect to an invariant world 



b ; time r [2]) 



and 



dK 1 

x v = = —9fivP (2) 

dK 1 dg^ v 

V --q — = ~~ 7; 3 P P 3 



lead to the geodesic equantion 

x p Tp Xi/Xp, (4) 
where what has appeared as a compatible connection form is given by 

1 ,dg>» dg*» dg»» 
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These results are tensor relations over the usual diffeomorphisms admitted by the 
manifold {x^}; writing the Hamiltonian in terms of (2), we see that the invariant interval 
on an orbit is proportional, through the constant Hamiltonian, to the square of the world 
time of evolution on the orbit, i.e., 

ds 2 = —Kdr 2 . (6) 
m 

We shall first study, in the following, a generalization of (1) consisting of the addition 
of a scalar field $>(x). The presence of such a scalar field can be considered as associated 
with the gauge covariant generalization of (1) in the Stuckelberg-Schrodinger equation [3] 
in the absence of four-vector gauge fields, an energy distribution not directly associated 
with visible light. We then show that there is a corresponding Hamiltonian K with a 
conformally modified metric, and no explicit additive scalar field, which has the form of 
Bekenstein's construction[4] for the realization of Milgrom's MOND program (modified 
Newtonian dynamics) [5] for achieving the observed galactic rotation curves. This simple 
form of Bekenstein's theory (called RAQUAL), which we discuss in detail in this work, for 
the sake of simplicity and clarity in the development of the mathematical method, does 
not properly account for causality and gravitational lensing; the theory has been further 
developed to include vector fields as well( TeVeS) [6] which has been relatively successful in 
accounting for these problems. We have shown previously that a gauge type Hamiltonian, 
with Minkowski metric and both vector and scalar fields [3] results, under a conformal 
map, in an effective Kaluza-Klein theory [7], and we shall indicate here (using a general 
Einsten metric) how the TeVeS structure can emerge, in terms of a Kaluza-Klein theory, 
in this way. More detailed analysis will be given in a subsequent publication. 

In the case treated in detail here, known as RAQUAL, the correspondence between 
K and K implies a relation between the conformal factor in K and the world scalar field 
and thus a possible connection between the so-called dark matter problem and a dark 
energy distribution represented by Application of the TeVeS theory can, furthermore, 
provide information on the Hamiltonian vector fields. 

2. Addition of a scalar potential 

The addition of a scalar potential to the Hamiltonian (1), in the form 

K=^-g^p» + <f>(x), (7) 

leads, according to the Hamilton equations, to the geodesic equation 1 

1 <9$ 

x p = — x v x^ — ~9p v ~Q^r- (8) 



Now, consider the Hamiltonian 



K = ^WV, (9) 



1 Note that (8) does not admit an equivalence principle. 
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where 

9nu = 4>gnv (10) 

This Hamiltonian can be put into correspondence with (7), as in the nonrelativistic case 
treated in [8], by defining 

</>=— ^- (11) 

with the constant mass shell constraint 

k = K = K. (12) 
As for (4), the Hamilton equations applied to (9) lead to the geodesic equation 2 

Xp = x u Xpj (13) 

where 

We remark that the construction based on Eqs. (9) and (10) admits the same family of 
diffeomorphisms as that of (7), since <p is scalar. Under these diffeomorphisms, both g pi/ 
and g^y are second rank tensors, and by construction of the connection forms, (4) and 
(13) are covariant relations. In the special coordinates for which (10) is taken explicitly, 
we have 

4— = ir-9 x '' + 4>4— , (15) 

ox v ox v ox v 

so that 



9 9 26 X dx v 9 d XfI 9 



(16) 



Substituting (11) into (16), this becomes, 

i , a$ <9$ 
9 p 2{k-$) x dx v p d Xfl 9 

and therefore the geodesic equation takes the form 

d$ 1 <9$ 



(17) 



x p - Vyx v x» m 9px dxx + k _^ dxv xpx v , (19) 



2 Eq.(13) does admit an equivalence principle, since g pu and T^ v are compatible. 
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This result (19) differs from the geodesic equation obtained from the Hamiltonian function 
K of Eq.(7). Let us, however, define a new velocity field, following the procedure used in 
[8], 

r = gT± v . (20) 

Solving for x v and substituting into the general form (13), with the identity 

we find the geodesic formula for the new velocity field 

r = -M&y V, (22) 

where 



- lit- < 23 > 



We assert that this result achieves a geometrical embedding of the motion generated 
by the Hamiltonian (7). Our method was to contruct the Hamiltonian (9) which generates 
geodesic equations with a compatible connection, thus providing a geometric basis for the 
theory. By the same methods used to test stability of orbits as used in ref. [8] , the geodesic 
deviation computed from the result (22) is effective in determining stability of the motion 
generated by the Hamiltonian (7). Applications of this type will be treated in a separate 
publication. To show that (22) is indeed a geometric embedding of the Hamiltonian (7), 
let us substitute the explicit form (10) for g^ v into (23). 

Using the definition (11), 

dg uX _ dcf) x dg v \ 



and the fact that 



one obtains 



dx^ dx^ dx^ 



= -T^Ti — y y — 0* — • 24 

2 oXfj, m ox^ 

Now, considering our transformation of velocity fields (20) heuristically as a local change 
of variables 3 , so that 

= s w A = <rV M #7, (25) 



dx^ dyP dyP ' 



3 The relation dy^ = gP v dx v is not integrable, and therefore does not uniquely define 
a set of coordinates For example, the "derivatives" = g v ^-^- are not commu- 

tative. The results of this identification have been, however, rigorously justified through a 
transformation of the affine parameter on the geodesic curves [9] . 
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we obtain 

r = -\g p »i^ry x --9> xp ^- p . (26) 

2 oyP m oy p 

This result differs in its structure from (8) in that the connection form contains just one 
term, while the connection form in (8) has three terms. To complete the equivalence, we 
define a related velocity field within the framework of the geodesic motions Let us 

define yet another velocity field 

z v = 9v^ (27) 

The derivative of y M then introduces an additional term, with indices symmetrized due to 
the bilinear form generated in the velocities. Using the relation (21) again (for derivatives 
of # M!y ), and identifying heuristically, in the same way as done above, 

d d 

dy~P =9p ^ (28) 

it follows from (26) that 

Zp = — T z^Za 9pv~^ ? (29) 
H m oz v 

where Y^ u is computed (in the same form) with all derivatives taken with respect to the 
variables {zp\ in place of the {xp\ in (5). Therefore, up to the transformation (27) within 
the family of velocity fields generated by the Hamilton equations from the conformally 
modified Hamiltonian (9), the geodesic equations (22) form a geometrical embedding of 
the original equations (8). Since T^ u and g^ v are compatible, there is a local flat space on 
this manifold in which parallel transport can be defined, and the tensor properties carry 
the same class of diffeomorphisms implicit in (7) and (8). 

We remark that the sequence of transformations(20) and (27) consists of 

z v = 9v^ = 9vp9^ X x\ = </> _1 iv, (30) 

independently of the coordinate system, since any Jacobians applied to these tensors will 
cancel. However, it is Eq. (22) that constitutes a nontrivial embedding of the orbits 
generated by (7). Our interest in this Section has been in relating the Hamiltonian (9) to 
the simplest Bekenstein-Milgrom form of MOND, without concern in the development of 
this simplified case for lensing or causal effects, for which a TeVeS type theory would be 
required, and with this, to be able to state restrictions on the form of the scalar field <E>. 
In the next Section, we indicate how a TeVeS can be generated in this framework, i.e., as 
a result of a conformal map. 

3. TeVeS and Kaluza-Klein Theory 

In this section, we show that the Te VeS theory can be cast into the form of a Kaluza- 
Klein construction. There has recently been a discussion[7], from the point of view of 
conformal correspondence, of a relativistic Hamiltonian with gauge invariant form 

K = - ea P){j> v - ea v ) -ea 5 , (31) 
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where the {a M }, as fields, may depend on the affine parameter r as well as x M , and the a 5 
field is necessary for the gauge invariance of the r derivative in the quantum mechanical 
Stueckelberg-Schrodinger equation. Here rj^ is the Minkowski metric (— 1, +1, +1, +1). 
As remarked in this work, Wesson[10] and Likofll] , as well as previous work on this struc- 
ture[3], have associated the a 5 field with mass density. It was shown[7] that a Hamiltonian 
of the form 

k=±-g^-ea»){ P »-ea») (32) 
can be put into correspondence with K by taking g^ u to have the conformal form 

k 

where k is the common (constant) value of K and K. In this correspondence, the equations 
of notion generated by K through the Hamilton equations, have extra terms beyond those 
provided by the connection form associated with g^, due to the presence of the gauge fields. 
Calculating the geodesic deviation, one could identify a curvature form associated with an 
effective five dimensional metric, consistent with the connection form in what then becomes 
the geodesic equation for the motion of a particle generated by the Hamilton equations 
obtained from K. This five dimensional effective metric is that of a Kaluza-Klein theory. 
We may apply the same procedure to the Hamiltonian 

K = ^9^ ~ WW - U v ) + <&, (34) 

where g^ is the Einstein metric, and $ is a world scalar field, and are gauge-like vector 
fields, as in Eq. (31). We shall give a more complete discussion of the dynamical properties 
of the equivalence in a subsequent paper, but it suffices for our purpose here to define, as 
in Eq. (33), the conformally modified metric 

~ ^ k 

9^-9^^—^ (35) 
= e ^ 9iiv'-, 

a Kaluza-Klein effective metric then emerges from the Hamilton equations applied to the 
"equivalent" Hamiltonian 

K=±-g^-WW -W), (36) 

as in ref.[7] 4 . 

Consider the Hamiltonian 

K K = ^J^V\ (37) 



4 One can choose _ftT(as in (32)) to be m/2, which results, according to the Hamilton 
equations, in dr as the invariant interval. 
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where [6] 

Qfjbv = e~ 2 Hg^ + U,Mu) - e 2< ^4 (38) 
This Hamiltonian then has the form 

K K = e- 2 %„p»p v - 2sinh2</>(W^) 2 , (39) 

where W M = g^JJ v ', i.e. with the same tensor properties as the fields appearing in Eq. 
(34). 

Let us now define a Kaluza-Klein type metric (of the form obtained in [7]) 

(40) 

Contraction to a bilinear form with the (5D) vectors pa = {p A ,P 5 }, with indices A = v on 
the right and A = \i on the left, one finds 

9abP A P B = g^p" + 2p 5 (p^) + (p 5 ) 2 g 55 . (41) 

If we take 

p 5 = -sinh20(^W M ) (42) 

and (755 = (the null choice of the constant assumed in ref.[7]), one sees that the Hamil- 
tonian (31) can be represented in terms of this Kaluza-Klein metric as 

K K = ^9ABP A P B - (43) 

Note that with the constraint that the fields U^ 1 are timelike unit vectors [6], enforced 
by using a Lagrange parameter, the product (p^U^) corresponds, in an appropriate local 
frame, to the energy of the particle, close to its mass in the case of a nonrelativistic particle, 
or to the frequency in the case of on-shell photons. It clearly remains to understand more 
deeply the apparently ad hoc choice of p 5 in (42) in terms of a 5D canonical dynamics, along 
with the structure of the 5D Einstein equations for g^B that follow from the geometry 
associated with (43). 



4. Conclusions 

A map of the type discussed in refs. [7], [8], of a Hamiltonian containing an Ein- 
stein metric, generating the connection form of general relativity, and a world scalar field, 
representing a distribution of energy on the spacetime manifold, into a corresponding 
Hamiltonian with a conformal metric (and compatible connection form), can account for 
the structure of the RAQUAL theory of Bekenstein and Milgrom[4]. Furthermore, ap- 
plying this correspondence to a Hamiltonian with gauge-type structure, we have shown 
that one obtains a non-compact Kaluza-Klein effective metric which can account for the 
TeVeS structure of Bekenstein, Sanders and Milgrom[6]. This method can be applied to 
the Brans-Dicke theory or other scalar-tensor theories as well. 
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The phenomenological constraints placed on the TeVeS variables in its astrophysical 
applications and by its MOND limit would, in principle, place constraints on the vector 
and scalar fields appearing in the corresponding Hamiltonian model. 
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Abstract 

It has been shown that the orbits of motion for a wide class of non-relativistic Hamil- 
tonian systems can be described as geodesic flows on a manifold and an associated dual by 



O 

(N 



means of a conformal map . This method can be applied to a four dimensional manifold 
of orbits in spacetime associated with a relativistic system. We show that a relativis- 
tic Hamiltonian which generates Einstein geodesies, with the addition of a world scalar 
field, can be put into correspondence in this way with another Hamiltonian with confor- 
Q-e mally modified metric. Such a construction could account for part of the requirements of 
Bekenstein for achieving the MOND theory of Milgrom in the post-Newtonian limit. The 
constraints on the MOND theory imposed by the galactic rotation curves, through this 
correspondence, would then imply constraints on the structure of the world scalar field. 
We then use the fact that a Hamiltonian with vector gauge fields results, through such a 
conformal map, in a Kaluza-Klein type theory, and indicate how the TeVeS structure of 
Bekenstein and Saunders can be put into this framework. We exhibit a class of infinitesimal 
! gauge transformations on the gauge fields U^(x) which preserve the Bekenstein- Sanders 
condition U^U 11 = — 1. The underlying quantum structure giving rise to these gauge fields 
is a Hilbert bundle, and the gauge transformations induce a non-commutative behavior to 
the fields, i.e. they become of Yang-Mills type. Working in the infinitesimal gauge neigh- 
ed borhood of the initial Abelian theory we show that in the Abelian limit the Yang-Mills 
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field equations provide residual nonlinear terms which may avoid the caustic singularity 
found by Contaldi et al. 

PACS 04.50.Cd,04.50.Kd,98.62.-g,95.35.+d,95.30.Sf 

1. Introduction 

The Hamiltonian [1] 

K = ^ V P^ (1) 

with Hamilton equations (written in terms of derivatives with respect to an invariant world 
time r [2]) 

op^ m 

and 

. _ dK 1 dg^ 



1 




(4) 



where what has appeared as a compatible connection form T p 



is given by 




1 p a/%m dg\v_ _ dg„ v , 
2 y { dx v dx» dx x) ' 



(5) 



These results are tensor relations over the usual diffeomorphisms admitted by the 
manifold writing the Hamiltonian in terms of (2), we see that the invariant interval 

on an orbit is proportional, through the constant Hamiltonian, to the square of the world 
time of evolution on the orbit, i.e., 



We shall first study, in the following, a generalization of (1) consisting of the addition 
of a scalar field $>(x). The presence of such a scalar field may be considered as a gauge 
compensation field for the r derivative in the evolution term of the covariant generalization 
of (1) in the Stuckelberg-Schrddinger equation [3], an energy distribution not directly 
associated with electromagnetic radiation in the usual sense. We then follow the method 
of ref. [4] to show that there is a corresponding Hamiltonian K with a conformally modified 
metric, and no explicit additive scalar field, which has the form of the construction of 
Bekenstein and Milgrom [5] for the realization of Milgrom's MOND program (modified 
Newtonian dynamics) [6] for achieving the observed galactic rotation curves. This simple 
form of Bekenstein's theory (called RAQUAL), which we discuss in some detail below 
for the sake of simplicity and clarity in the development of the mathematical method, 
does not properly account for causality and gravitational lensing; the theory has been 
further developed to include vector fields (which we shall call Bekenstein-Sanders fields) as 
well(TeVeS) [7], which has been relatively successful in accounting for these problems. It 
has been shown [8], moreover, that a gauge type Hamiltonian, with Minkowski metric and 
both vector and scalar fields results, under a conformal map, in an effective Kaluza-Klein 
theory. We shall indicate here (using a general Einstein metric) how the Te VeS structure 
can emerge in terms of a Kaluza-Klein theory in this way. It is essential in this construction 
that the Bekenstein-Sanders fields be considered as gauge fields. As a realization of this 
possibility, we exhibit an infinitesimal gauge transformation on the underlying quantum 
theory for which the vector fields, which we shall call U^{x), emerge as gauge compensation 
fields, such that, as required by the TeVeS theory the property U^U^ = — 1 [7] is preserved 
under such gauge transformations. The corresponding quantum theory then has the form 
of a Hilbert bundle and, in this framework, the gauge fields are of (generalized) Yang-Mills 
type [9]. Working in the infinitesimal neighborhood of a gauge in which the fields are 
Abelian, we show that in the limit the contributions from the nonabelian sector provide 
nonlinear terms in the field equations which may avoid the caustic singularity found by 
Contaldi et al [10]. Further investigation of this structure will be given in a subsequent 
publication. 

For both the RAQUAL and the TeVeS theories, the correspondence between K and 
K implies a relation between the conformal factor in K and the world scalar field and 



ds 2 = —Kdr 2 . 



(6) 



rn 
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thus a possible connection between the so-called dark matter problem and a dark energy 
distribution represented by <E>. 

2. Addition of a scalar potential and conformal equivalence 

The addition of a scalar potential to the Hamiltonian (1), in the form 

K=^-g^ Pfl p u + ^(x), (7) 

leads, according to the Hamilton equations, to the geodesic equation 1 

1 <9$ 

& = -T% v x»x» - -g pv j—. (8) 
m ox" 

Now, consider the Hamiltonian (we carry out the calculations explicitly here since we 
shall have need of some of the intermediate results) 

K = 7^9^(y)p»P». (9) 
It follows from the Hamilton equations that 



OK 



so that 

Vv = mg^y* 1 ( 10 ) 

and 

_dK 1 dg^ 

As in (4), it then follows that 

r = -Ky x y a , (ii) 

where, as for (4), 

f*i _ 1 ^ f d 9™ i d<jv\ ggAa -i n9 v 
Act " 2 9 1 «V <%/ CT c>r 1 J 

We now establish an equivalence between the Hamiltonians (7) and (9) by assuming 
the momenta p^ equal at every moment r in the two descriptions. With the constraint 

K = K = k, (13) 

if we assume the conformal form 

g v °{y) = <P(y)g^(x), (14) 



1 Note that (8) does not admit an equivalence principle, but (11), arising from (9) does. 

3 



it follows that 

4(y)(k-$(x)) = k. (15) 

The relation (15) is not sufficient to construct y as a function of x, but if we impose 
the relation 

Sx" = ( f>- 1 (y)8y t * (16) 

between variations generated on position in the two coordinate systems, it is sufficient to 
evaluate derivatives of <p(y) in terms of derivatives with respect to x of the scalar field 
$(x)[ll](see also [12]). We review this construction briefly below. 

We remark that the construction based on Eqs. (9) and (14) admits the same family 
of diffeomorphisms as that of (7), since (p is scalar. Under these diffeomorphisms, both 
g^ v and are second rank tensors, and by construction of the connection forms, (4) and 
(11) are covariant relations. 

To see how these derivatives are constructed on the constraint hypersurface determined 
by (15), let us, for brevity, define 

f <*> s *tW (l7) 

so that the constraint relation (15) reads 

<t>(y) = F(x). (18) 
Then, since variations in x and y are related by (16), 

dF(x) 

(j)(y + 5y) = F(x + 5x)^F(x) + 5x^-^. (19) 



To first order in Taylor's series on the left, we obtain the relation 

(20) 



dy^ dx^ 
We may therefore define a derivative, restricted to the constraint hypersurface 

s _m = ^ {y) ?m (21) 

dy» r KyJ dx» V ' 

The Leibniz relation follows easily for the product of functions, it e.g., for (p(y)g^ u (x). 

In a similar way, the second derivative can be obtained from (19) by recognizing that 
the variation in x is to be computed at the point y + Sy. Keeping terms of second order in 
the expansion of both sides, one can define the second derivative restricted to the constraint 
hypersurface defined by (18); although it appears that a second derivative defined by (21) 
would not be symmetric, both the derivative of (21) and the second derivative computed 
from (19) on the constraint hyperfurface agree and are symmetric [11], i.e., 

~d*F(x) = d 2 F(x) 
dy^dy v dy v dy^ 
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This implies that the restricted derivative defined by (21) behaves as a bona fide derivative 
on the constraint hypersurface, admitting the consistent coexistence of the coordinates x 
and y related by (15). We will not have further use of (22) here, primarily relevant for the 
calculation of stability criteria through geodesic deviation (the second derivative occurs in 
the curvature tensor). 

In the following, we complete our argument of equivalence by reconstructing the equa- 
tions of motion following from the Hamilton equations applied to (7), i.e., Eq. (8). 

We begin our reconstruction, in analogy with the procedure used in the nonrelativistic 
problem [8], by defining the new variable such that 

z» = 9pu(y)y u - (23) 

Substituting y u = g pu (y)z^ into (11), the r derivatives of g pu {y) generate terms that cancel 
two of the terms in T^ a , leaving 

Z »=2W VV ' ( ) 

Now, substituting for y x from (23), and using the identity 

dg\a . ap = dg 
dy u dy 



9*%9" = -??. (25) 



we find 

1 d(p p 

* v = ~2~dr iyip ' (26) 

Finally, from the variational type argument we used above, 



gP^ y + Sy)-g^(y) = -f^Sy" 



dy» 
so that 



pi 

g vX 6z> 



or 





QgPl 




dz x 


QgPl ^ 


QgPl 




dz x 



We therefore have the alternative form 



1 dg p ^ 

z» = ~2 9uX ~dzx~ ZpZl ' ( 29 > 



(27) 



(28) 
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This result constitutes a "geometric" embedding of the Hamiltonian motion induced by 
(7) in the same way as for the nonrelativistic case. Substituting the explicit form of g pl 
in terms of the original Einstein metric from (14), one obtains 

1 dg^. . 1 x dcf> . 
z v = ~2 9uX ~dz^ zpZl ~ 2 9vX ~dz~ x g lZp ' ' 

The second term contains the potential field, as in the Hamilton equations, but the first 
term does not contain the full connection form. We may finally, however, define a "decon- 
traction" of the connection in (30) using the Einstein metric. In fact, since according to 
(16), r = and by (23), 

z» = g^y u = <p~ 1 g ll vy u , (31) 

it follows that 

&n = 9^x v . (32) 

Making this substitution in (30) leads explicitly, taking into account the k shell constraint 
(13) and the form of (7), to the equation (8). We have thus completed our demonstration of 
the equivalence between the purely metric form of the Hamiltonian (9) and the Hamilton 
(7), for which the relation (29) corresponds to a dynamics generated by a compatible 
connection form, and constitute a "geometric" embedding of the original Hamiltonian 
motion. 

Our interest in this section has been in relating the Hamiltonian (7) to the simplest 
Bekenstein-Milgrom form of MOND, without concern in the development of this simplified 
case for lensing or causal effects, for which a TeVeS type theory would be required. In the 
next Section, we indicate how a TeVeS theory can be generated in this framework, i.e., 
as a result of a conformal map. 



3. TeVeS and Kaluza-Klein Theory 

In this section, we show that the Te VeS theory can be cast into the form of a Kaluza- 
Klein construction. There has recently been a discussion [8], from the point of view of 
conformal correspondence, of the equivalence of a relativistic Hamiltonian with an electro- 
magnetic type gauge invariant form [3] (here 7f~ v is the Minkowski metric (—1, +1, +1, +1)) 

K = -^V^iPti ~ ea^)(p u - ea u ) - ea 5 , (33) 

where the {a M }, as fields, may depend on the affine parameter r as well as x M , and the a 5 
field is necessary for the gauge invariance of the r derivative in the quantum mechanical 
Stueckelberg-Schrodinger equation, with a Kaluza-Klein theory. As remarked in this work, 
Wesson [13] and Liko [14] , as well as previous work on this structure[3], have associated 
the source of the 05 field with mass density. A Hamiltonian of the form 

K = ^-g* 1 "^ - eaj(p v - ea u ) (34) 
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can be put into correspondence with K by taking g^ u to be 

r = ^ * (35) 
k + eas 

where k is the common (constant) value of K and K. In this correspondence, the equations 
of notion generated by K through the Hamilton equations, have extra terms, beyond those 
provided by the connection form associated with g^ v ', due to the presence of the gauge fields. 
These additional terms can be identified as belonging to a connection form associated with 
a five dimensional metric, that of a Kaluza-Klein theory. 
We may apply the same procedure to the Hamiltonian 

K = ^9^(P» ~ <MJ(P» ~ <Mv) + ^ (36) 

where g^ v is an Einstein metric, $ is a world scalar field, and W M are to be identified with 
the Bekenstein- Sanders fields for which[6] U U U U = — 1, with = g^ v U v . 

We discuss in Section 4 a class of gauge transformations on the wave functions of the 
underlying quantum theory for which the arise as gauge compensation fields. 

Let us define, as in Eq. (35), the conformally modified metric 



k-$ (37) 

= e" 5 

The "equivalent" Hamiltonian 



K = ^""(Pv ~ <Mj(Pv ~ <Mv) (38) 

then generates, through the Hamilton equations, an equation of motion which corresponds 
to the geodesic equation for an effective Kaluza-Klein metric, as in ref.[8]. 
Now, consider the Hamiltonian 

K K = ^'"'PvPu, (39) 

with the Bekenstein- Sanders metric[7] 

~ gi xv = e -24>( g »» + _ e i<t>um v (40) 

The Hamiltonian Kk then has the form 

K K = e-^g^p^ - 2 sinh 2</>(W> M ) 2 , (41) 

Let us now define a Kaluza-Klein type metric of the form obtained in [7] , arising from 
the equations of motion generated by (38), 

«"-(£ W- (42) 
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Contraction to a bilinear form with the (5D) vectors pa = {px^Ps}, with indices A = v on 
the right and A = \x on the left, one finds 



9 Ab papb = g^P^Pu + 2 P5 (p^) + ( P5 ) 2 g 



(43) 



If we take 




(1 ± y/l - 2# 55 sinh2(/>) , 



(44) 



then the Kaluza-Klein theory coincides with (41), i.e., 



K K = - — 9 PaPb- 



(45) 



As discussed by Wesson [13], Kaluza [15] chose (755 = const, for consistency with electro- 
magnetism, while Wesson [13] makes the more general choice of a world scalar field. In 
particular, the value g 55 = is well defined (as in [8]). 

Since the fields are timelike unit vectors [7], (p^U^) corresponds, in an appropriate 
local frame, to the energy of the particle, close to its mass in the case of a nonrelativistic 
particle, or to the frequency in the case of on-shell photons. It clearly remains to understand 
more deeply the apparently ad hoc choice of p 5 in (44) in terms of a 5D canonical dynamics, 
along with the structure of the 5D Einstein equations for qab that follow from the geometry 
associated with (45). We shall study these questions in a succeeding paper. 

4 The Bekenstein- Sanders Vector Field as a Gauge Field 

Essential features of the Bekenstein- Sanders field [7] of the TeVeS theory are that it 
be a local field, i.e., U^ix), and there is a normalization constraint 



so that the vector is timelike. To preserve the normalization condition (46) under gauge 
transformation, we shall study the construction of a class of gauge transformations which 
essentially moves the U{x) field on a hyperbola with a Lorentz transformation (at the point 



If we think of our underlying quantum structure, which generates the gauge field, as 
a fiber bundle with base x M , then we must think of the transformation acting in such a 
way that the absolute square (norm) of the wave function attached to the base point x M 
preserves its value [9]. 

An analogy can be drawn to the usual Yang-Mills gauge [9] on SU(2), where there 
is a two- valued index for the wave function ip a (x). The gauge transformation in this case 
is a two by two matrix function of x, and acts only on the indices a. The condition of 
invariant absolute square (probability) is 



WU» = -1, 



(46) 



x). 




(47) 



8 



Generalizing this structure, one can take the indices a to be continuous, so that (47) 
becomes 

J{dU)\ j \dU')U{UM')^{U\x)\ 2 = j \dU)\i){U,x)\ 2 , (48) 

implying that U(U,W) is a unitary operator on a Hilbert space L 2 (dU). Since we are 
assuming that lies on an orbit determined by (48), the measure is 

= (49) 

i.e., a three dimensional Lorentz invariant integration measure. 

Moreover, the Lorentz transformation on is generated by a non-commutative oper- 
ator, and therefore the gauge transformation is non-Abelian. We demonstrate the resulting 
noncommutativity of the operator valued fields, W, after an infinitesimal gauge transfor- 
mation of ths type, explicitly below. 

This construction is somewhat similar to the treatment of the electromagnetic poten- 
tial vector and its time derivative as oscillator variables in the process of second quantiza- 
tion of the radiation field (the energy density of the field is given by these variables in the 
form of an oscillator). 

We now examine the gauge condition: 

( Pfl -eU^U^ = U( Pfl -eU^ (50) 
Identifying Pil with —id/dx^, and cancelling the terms Up^ifj on both sides, we obtain 

W, = UU,V-i - ~U~\ (51) 

in the same form as the Yang-Mills theory [9] . It is evident in the Yang-Mills theory, that 
due to the matrix nature of the second term, the field will be algebra-valued, resulting in 
the usual structure of the Yang-Mills non-Abelian gauge theory. Here, if the transformation 
U is a Lorentz transformation, the numerical valued field would be carried, in the first 
term, to a new value on a hyperbola. However, the second term may well be operator 
valued on L 2 (dU), and thus, as in the Yang-Mills theory, W' M would become nonabelian. 
It follows from (51) that the field strengths 

! ^- 9 ^- S ^ + w -^ ] (52) 

are related to the the field strengths in the transformed form 



according to 

fU x ) = Uf^{x)U~\ (54) 
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just as in the finite dimensional Yang-Mills theories. 
This result follows from writing out, from (51), 

dU' 8U , dU u , OU' 1 

i d 2 U u _ 1 i dU dU- 1 
~ e dx^dx v ~ edx^ dx u ' 

and subtracting the same expression with /i, v reversed. Then add the result to 



QTJ 

r dU , i r dU i dU 

Whenever the combination 
appears, it should be replaced by 

ou- 1 



(56) 



<9a^ 

The result (54) then follows after a little manipulation. 

Now, consider the possibility that this finite gauge transformation leaves U^U^ = — 1. 
We write out 

(lIU^U-i - ^U-^UWU^ - l -^U~i) = -1 - igLu»U-i 
e ox^ e ax^ e ax^ 

i , dU , 

= -i _ « + ^ af/_1 



e cte^ e M dx^ 

^_3U_dU^_ 
e 2 <9a; M <9x M 

(57) 

It may be possible that U can be chosen to make all but the first term in (57) vanish, 
but in the case of finite gauge transformations, it is not so easy to see how to construct 
examples. For the infinitesimal case, it is, however, easy to construct a gauge function 
with the required properties. For 

U = l + iG, (58) 

where G is infinitesimal, (51) becomes 

1 dG 



U^=U,+ i[GM»] + + 0(G 2 ). (59) 
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Then, 



Let us take 



c = -?EKw*).^-^)} 



2 ^ L <9W 7 <9W A / 61 \ 



= -X){wA 7 (W,x),^} 



2 

where symmetrization is required since is a function of U as well as x, and 

This construction is valid in the initially special gauge, which we shall call the "special 
abelian gauge", in which the components of W M commute. The appearance of VC 1 in 
the gauge functions is then admissible since this quantity acts on the wave functions 
< U,x\ip) = ip{U,x) at the point x,in the representation in which the operator W M on 
L 2 (dU) is diagonal. 

Our investigation in the following will be concerned with a study of the infinitesi- 
mal gauge neighborhood of this limit, where the components of W M do not commute, and 
therefore constutite a Yang Mills type field. We shall show in the limit that the corre- 
sponding field equations acquire nonlinear terms, and may therefore suppress the caustic 
singularities found by Contaldi et al [10]. They found that nonlinear terms associated with 
a non-Maxwellian type action, such as (<9 M ZY^) , could avoid this caustic singularity, so 
that the nonlinear terms we find as a residue of the Yang-Mills structure induced by our 
gauge transformation might achieve this effect in a natural way. 

The second term of (60), which is the commutator of G with U^U^ vanishes, since 
this product is Lorentz invariant (the symmetrization in G does not affect this result). 

We now consider the third term in (60). 



= \\N x "< dUXl W + dUXl iV A W (63) 
2 L dx^ dx^ 

dx* dx» J 



There are two terms proportional to 



dx» 

If we take (locally) 



wa 7 (W, x) = u\y(k u x u ), (64) 
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where k v U v = 0, then 

^-U, = k.WJ Xl = 0. (65) 

For the remaining two terms, 



dx^ dx^ 



= N X ^W 



,du Xl du Xl . 



(66) 



+ [U», N^] + zr^LWN* 
1 J dx^ dx» 

Since the commutators contain only terms linear in W M and they have opposite sign, they 
cancel. The remaining terms are zero by the argument (65). The condition U^W 1 = — 1 is 
therefore invariant under this gauge transformation, involving the coefficient u Xl which is a 
function of the projection of rr M onto a hyperplane orthogonal to U^, i.e., a function of 
where kyJA^ = 0. The vector of course, depends on^ (for example, k^ = 14^(14 •b) + b /J ,, 
for some 6 M 7^ 0). 

We now demonstrate explicitly the nonabelian nature of the gauge fields after in- 
finitesinal gauge transformation. With (59), the commutator term in (53) is 

1 3C 1 3C 

MM = (W„ + *[G,W M ] + -q^Mu + i[G,U v ] + 

1 c)C 1 r)C 

- (M „ + i[G , w j + -_)( M , + i[G , w j + -_) (67) 

+ z[W M ,[G,W„]]-z[W,,[G,W M ]], 

where the remaining terms have identically cancelled. Note that this expression does not 
contain any noncommutative quantities. Now, 

[G, U v \ = 2ieu^Uj (68) 

and 



The terms involving [G, U v \ and [G, ZY M ] therefore cancel, so that 

dx v dx^ 



KM = 2%(^-tf) (™) 



We have taken u) X fl = uj x ^(k a x a ), so that 

du xf " 



dx v 
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= (71) 



and therefore 

[Ul,U' v ] = 2i(Ku'\ - k^'\)U x , (72) 

generally not zero. This demonstrates the nonabelian character of the fields. In the Abelian 
limit, we may take u' — > 0, but as we shall a residual nonlinearity, which depends on u" 
may remain in the field equations 

We now consider the derivation of field equations from a Lagrangian constructed with 
the ip's and f^ v f^ v . We take the Lagrangian to be of the form (the indices are raised and 
lowered with g^ v ) 

C = C f + C m , (73) 

where 

£f = -\rf»» (74) 

and 

Q 1 

C m = r(ifr-^(Pv-eU l *)9' tv (Pu-eU v )-*)l> + c.c. (75) 

We shall be working in the infinitesimal neighborhood of the special gauge for Abelian 
Ufx, for which it has the form given in (59) for infinitesinal G. It is therefore not Abelian to 
first order, but we take its variation 8U to be a c-number function, carrying the variation, 
to lowest order, by variation of the first term in (59), and not varying the part of U 
introduced by the infinitesimal gauge transformation (evaluated on the original value of 

u). 

In carrying out the variation of C m , the contributions of varying the ^'s with respect 
to U vanish due to the field equations (Stueckelberg-Schrodinger equation) obtained by 
varying ift* (or ift), and therefore in the variaton with respect to U, only the explicit 
presence of U in (75) need be taken into account. 

Note that for the general case of U generally operator valued, we can write 

r(p» - dA^{p u - tU v )i> = g^((p» - eU^YiPu - d4 v )4>, (76) 

since the Lagrangian density (75) contains an integration over (dU')(dU") (considered in 
lowest order) as well as an integration over (dx) in the action and the operators U are 
Hermitian. In the limit in which U is evaluated in the special Abelian gauge (real valued) , 
and noting that p M is represented by an imaginary differential operator, we can write this 

as 

^VCpa. - MJipv - = -g^ip* + <U»Wip v - eU„)ii, (77) 

i.e., replacing explicitly p^ by —i^d/dx^) = — i<9 M , we have 

SuC m = -^{V>*(<9„ - ieU»)if> - ((^ + ieU^r)^}SU^, (78) 
where we have called g^ v bU v = 5U^, or, 

S u C m = j^U,x)SU^ (79) 
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where j M has the usual form of a gauge invariant current. 

For the calculation of the variation of Cf we note that the commutator term in (52) 
is, in lowest order, a c-number function, as given in (72). 

Calling 

u'\U x EE Vp, (80) 

we compute the variation of 



Then, for 



we compute 



[U' li ,U , v ]=2i(k v v lt -k li v v ) (81) 



SuKM = Su^KMl (82) 

KM = n^v, + K^f) - u)). (83) 



With our choice of k v = U V {U • b) + b v 

dk 



. )u - = 8S(U-b)+U„b\ (84) 



so that 



[W;, U' v \ = 2i{5^{U ■ b) + Uub^v" 



aw 7 



= O 1 



z.e. 



S u [U , li M = Cr iiV SU r (86) 

The quantity u p is proportional to the derivative of ui*. In the limit that ui,ui' — > 
(cf. (81)), the second derivative, u" which appears in O 1 ^ u may not vanish (somewhat 
analogous to the case in gravitional theory when the connection form vanishes but the 
curvature does not), so that this term can contribute in limit to the special Abelian gauge. 
Returning to the variation of Cf in (74), we see that 

6C f = -d v f^bW + 2if^S[U^,U v }, (87) 

where we have taken into account the fact that \U/j,,U u ] is a commuting function, and 
integrated by parts the derivatives of 5U. With (86) we obtain 

Act 



5C f = -d v j^W + 2ief XlT A \8U» 

Since the coefficient of SW 1 must vanish, we obtain, with (79), the Yang-Mills equations 
for the fields given the source currents 

d u U„ = j»-2ief Xa O x %, (89) 
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which is nonlinear in the fields U^, as we have seen, even in the Abelian limit, where, from 
(78) and (79), 

3„ = - idA^ - + idA^W)^}. (90) 

We point out that this current corresponds to a flow of the matter field; the absolute 
square of the wave functions corresponds to an event density. The coupling e is not 
necessarily charge, and the fields U are not necessarily electromagnetic even in the Abelian 
limit. However, the Hamiltonian (36) leads directly to a Lorentz type force, similar in 
form to that generated by the Hilbert-Einstein action. The dynamics of this system will 
be investigated in a forthcoming paper. 

5. Conclusions 

A map of the type discussed in ref. [8] , of a Hamiltonian containing an Einstein metric, 
generating the connection form of general relativity, and a world scalar field, representing 
a distribution of energy on the spacetime manifold, into a corresponding Hamiltonian with 
a conformal metric (and compatible connection form), can account for the structure of 
the RAQUAL theory of Bekenstein and Milgrom[5]. Furthermore, applying this corre- 
spondence to a Hamiltonian with gauge-type structure, we have shown that one obtains a 
non-compact Kaluza-Klein effective metric which can account for the TeVeS structure of 
Bekenstein, Sanders and Milgrom[7]. 

In order to maintain the constraint condition UpJA^ = — 1 for the Bekenstein-Sanders 
fields, under local gauge transformations, we have introduced a class of gauge of gauge 
transformations on the underlying quantum theory which acts on the Hilbert bundle, 
quite analogous to that arising in the second quantization of the electromagnetic field 
(where the vector potentials and their time derivatives are considered as quantum oscillator 
variables) associated with the values of the gauge fields. The action of this class of gauges 
induces a nonabelian structure on the fields, which therefore satisfy Yang-Mills type field 
equations with source currents associated with matter flow. In the Abelian limit, these 
equations contain residual non-linear terms which may avoid the caustic singularities found 
by Contaldi et al for an electromagnetic type gauge field. 

The phenomenological constraints placed on the TeVeS variables in its astrophysical 
applications and on its MOND limit [16] would, in principle, place constraints on the vector 
and scalar fields appearing in the corresponding Hamiltonian model, for which the additive 
world scalar field corresponds to an energy distribution not associated with electromagnetic 
radiation. 
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Abstract 

It has been shown that the orbits of motion for a wide class of non-relativistic Hamil- 
tonian systems can be described as geodesic flows on a manifold and an associated dual by 



O 

(N 



means of a conformal map . This method can be applied to a four dimensional manifold 
of orbits in spacetime associated with a relativistic system. We show that a relativis- 
tic Hamiltonian which generates Einstein geodesies, with the addition of a world scalar 
field, can be put into correspondence in this way with another Hamiltonian with confor- 
Q-e mally modified metric. Such a construction could account for part of the requirements of 
Bekenstein for achieving the MOND theory of Milgrom in the post-Newtonian limit. The 
constraints on the MOND theory imposed by the galactic rotation curves, through this 
correspondence, would then imply constraints on the structure of the world scalar field. 
We then use the fact that a Hamiltonian with vector gauge fields results, through such a 
conformal map, in a Kaluza-Klein type theory, and indicate how the TeVeS structure of 
Bekenstein and Saunders can be put into this framework. We exhibit a class of infinitesimal 
! gauge transformations on the gauge fields U^(x) which preserve the Bekenstein- Sanders 
condition U^U 11 = — 1. The underlying quantum structure giving rise to these gauge fields 
is a Hilbert bundle, and the gauge transformations induce a non-commutative behavior to 
the fields, i.e. they become of Yang-Mills type. Working in the infinitesimal gauge neigh- 
ed borhood of the initial Abelian theory we show that in the Abelian limit the Yang-Mills 



O 



field equations provide residual nonlinear terms which may avoid the caustic singularity 
found by Contaldi et al. 

PACS 04.50.Cd,04.50.Kd,98.62.-g,95.35.+d,95.30.Sf 

1. Introduction 

The Hamiltonian [1] 

K = ^ V P^ (1) 

with Hamilton equations (written in terms of derivatives with respect to an invariant world 
time r [2]) 

op^ m 

and 

. _ dK 1 dg^ 



1 




(4) 



where what has appeared as a compatible connection form T p 



is given by 




1 p a/%m dg\v_ _ dg„ v , 
2 y { dx v dx» dx x) ' 



(5) 



These results are tensor relations over the usual diffeomorphisms admitted by the 
manifold writing the Hamiltonian in terms of (2), we see that the invariant interval 

on an orbit is proportional, through the constant Hamiltonian, to the square of the world 
time of evolution on the orbit, i.e., 



We shall first study, in the following, a generalization of (1) consisting of the addition 
of a scalar field $>(x). The presence of such a scalar field may be considered as a gauge 
compensation field for the r derivative in the evolution term of the covariant generalization 
of (1) in the Stuckelberg-Schrddinger equation [3], an energy distribution not directly 
associated with electromagnetic radiation in the usual sense. We then follow the method 
of ref. [4] to show that there is a corresponding Hamiltonian K with a conformally modified 
metric, and no explicit additive scalar field, which has the form of the construction of 
Bekenstein and Milgrom [5] for the realization of Milgrom's MOND program (modified 
Newtonian dynamics) [6] for achieving the observed galactic rotation curves. This simple 
form of Bekenstein's theory (called RAQUAL), which we discuss in some detail below 
for the sake of simplicity and clarity in the development of the mathematical method, 
does not properly account for causality and gravitational lensing; the theory has been 
further developed to include vector fields (which we shall call Bekenstein-Sanders fields) as 
well(TeVeS) [7], which has been relatively successful in accounting for these problems. It 
has been shown [8], moreover, that a gauge type Hamiltonian, with Minkowski metric and 
both vector and scalar fields results, under a conformal map, in an effective Kaluza-Klein 
theory. We shall indicate here (using a general Einstein metric) how the Te VeS structure 
can emerge in terms of a Kaluza-Klein theory in this way. It is essential in this construction 
that the Bekenstein-Sanders fields be considered as gauge fields. As a realization of this 
possibility, we exhibit an infinitesimal gauge transformation on the underlying quantum 
theory for which the vector fields, which we shall call U^{x), emerge as gauge compensation 
fields, such that, as required by the TeVeS theory the property U^U^ = — 1 [7] is preserved 
under such gauge transformations. The corresponding quantum theory then has the form 
of a Hilbert bundle and, in this framework, the gauge fields are of (generalized) Yang-Mills 
type [9]. Working in the infinitesimal neighborhood of a gauge in which the fields are 
Abelian, we show that in the limit the contributions from the nonabelian sector provide 
nonlinear terms in the field equations which may avoid the caustic singularity found by 
Contaldi et al [10]. Further investigation of this structure will be given in a subsequent 
publication. 

For both the RAQUAL and the TeVeS theories, the correspondence between K and 
K implies a relation between the conformal factor in K and the world scalar field and 



ds 2 = —Kdr 2 . 



(6) 



rn 
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thus a possible connection between the so-called dark matter problem and a dark energy 
distribution represented by <E>. 

2. Addition of a scalar potential and conformal equivalence 

The addition of a scalar potential to the Hamiltonian (1), in the form 

K=^-g^ Pfl p u + ^(x), (7) 

leads, according to the Hamilton equations, to the geodesic equation 1 

1 <9$ 

& = -T% v x»x» - -g pv j—. (8) 
m ox" 

Now, consider the Hamiltonian (we carry out the calculations explicitly here since we 
shall have need of some of the intermediate results) 

K = 7^9^(y)p»P». (9) 
It follows from the Hamilton equations that 



OK 



so that 

Vv = mg^y* 1 ( 10 ) 

and 

_dK 1 dg^ 

As in (4), it then follows that 

r = -Ky x y a , (ii) 

where, as for (4), 

f*i _ 1 ^ f d 9™ i d<jv\ ggAa -i n9 v 
Act " 2 9 1 «V <%/ CT c>r 1 J 

We now establish an equivalence between the Hamiltonians (7) and (9) by assuming 
the momenta p^ equal at every moment r in the two descriptions. With the constraint 

K = K = k, (13) 

if we assume the conformal form 

g v °{y) = <P(y)g^(x), (14) 



1 Note that (8) does not admit an equivalence principle, but (11), arising from (9) does. 
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it follows that 

4(y)(k-$(x)) = k. (15) 

The relation (15) is not sufficient to construct y as a function of x, but if we impose 
the relation 

Sx" = ( f>- 1 (y)8y t * (16) 

between variations generated on position in the two coordinate systems, it is sufficient to 
evaluate derivatives of <p(y) in terms of derivatives with respect to x of the scalar field 
$(x)[ll](see also [12]). We review this construction briefly below. 

We remark that the construction based on Eqs. (9) and (14) admits the same family 
of diffeomorphisms as that of (7), since (p is scalar. Under these diffeomorphisms, both 
g^ v and are second rank tensors, and by construction of the connection forms, (4) and 
(11) are covariant relations. 

To see how these derivatives are constructed on the constraint hypersurface determined 
by (15), let us, for brevity, define 

f <*> s *tW (l7) 

so that the constraint relation (15) reads 

<t>(y) = F(x). (18) 
Then, since variations in x and y are related by (16), 

dF(x) 

(j)(y + 5y) = F(x + 5x)^F(x) + 5x^-^. (19) 



To first order in Taylor's series on the left, we obtain the relation 

(20) 



dy^ dx^ 
We may therefore define a derivative, restricted to the constraint hypersurface 

s _m = ^ {y) ?m (21) 

dy» r KyJ dx» V ' 

The Leibniz relation follows easily for the product of functions, it e.g., for (p(y)g^ u (x). 

In a similar way, the second derivative can be obtained from (19) by recognizing that 
the variation in x is to be computed at the point y + Sy. Keeping terms of second order in 
the expansion of both sides, one can define the second derivative restricted to the constraint 
hypersurface defined by (18); although it appears that a second derivative defined by (21) 
would not be symmetric, both the derivative of (21) and the second derivative computed 
from (19) on the constraint hyperfurface agree and are symmetric [11], i.e., 

~d*F(x) = d 2 F(x) 
dy^dy v dy v dy^ 
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This implies that the restricted derivative defined by (21) behaves as a bona fide derivative 
on the constraint hypersurface, admitting the consistent coexistence of the coordinates x 
and y related by (15). We will not have further use of (22) here, primarily relevant for the 
calculation of stability criteria through geodesic deviation (the second derivative occurs in 
the curvature tensor). 

In the following, we complete our argument of equivalence by reconstructing the equa- 
tions of motion following from the Hamilton equations applied to (7), i.e., Eq. (8). 

We begin our reconstruction, in analogy with the procedure used in the nonrelativistic 
problem [8], by defining the new variable such that 

z» = 9pu(y)y u - (23) 

Substituting y u = g pu (y)z^ into (11), the r derivatives of g pu {y) generate terms that cancel 
two of the terms in T^ a , leaving 

Z »=2W VV ' ( ) 

Now, substituting for y x from (23), and using the identity 

dg\a . ap = dg 
dy u dy 



9*%9" = -??. (25) 



we find 

1 d(p p 

* v = ~2~dr iyip ' (26) 

Finally, from the variational type argument we used above, 



gP^ y + Sy)-g^(y) = -f^Sy" 



dy» 
so that 



pi 

g vX 6z> 



or 





QgPl 




dz x 


QgPl ^ 


QgPl 




dz x 



We therefore have the alternative form 



1 dg p ^ 

z» = ~2 9uX ~dzx~ ZpZl ' ( 29 > 



(27) 



(28) 
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This result constitutes a "geometric" embedding of the Hamiltonian motion induced by 
(7) in the same way as for the nonrelativistic case. Substituting the explicit form of g pl 
in terms of the original Einstein metric from (14), one obtains 

1 dg^. . 1 x dcf> . 
z v = ~2 9uX ~dz^ zpZl ~ 2 9vX ~dz~ x g lZp ' ' 

The second term contains the potential field, as in the Hamilton equations, but the first 
term does not contain the full connection form. We may finally, however, define a "decon- 
traction" of the connection in (30) using the Einstein metric. In fact, since according to 
(16), r = and by (23), 

z» = g^y u = <p~ 1 g ll vy u , (31) 

it follows that 

&n = 9^x v . (32) 

Making this substitution in (30) leads explicitly, taking into account the k shell constraint 
(13) and the form of (7), to the equation (8). We have thus completed our demonstration of 
the equivalence between the purely metric form of the Hamiltonian (9) and the Hamilton 
(7), for which the relation (29) corresponds to a dynamics generated by a compatible 
connection form, and constitute a "geometric" embedding of the original Hamiltonian 
motion. 

Our interest in this section has been in relating the Hamiltonian (7) to the simplest 
Bekenstein-Milgrom form of MOND, without concern in the development of this simplified 
case for lensing or causal effects, for which a TeVeS type theory would be required. In the 
next Section, we indicate how a TeVeS theory can be generated in this framework, i.e., 
as a result of a conformal map. 



3. TeVeS and Kaluza-Klein Theory 

In this section, we show that the Te VeS theory can be cast into the form of a Kaluza- 
Klein construction. There has recently been a discussion [8], from the point of view of 
conformal correspondence, of the equivalence of a relativistic Hamiltonian with an electro- 
magnetic type gauge invariant form [3] (here 7f~ v is the Minkowski metric (—1, +1, +1, +1)) 

K = -^V^iPti ~ ea^)(p u - ea u ) - ea 5 , (33) 

where the {a M }, as fields, may depend on the affine parameter r as well as x M , and the a 5 
field is necessary for the gauge invariance of the r derivative in the quantum mechanical 
Stueckelberg-Schrodinger equation, with a Kaluza-Klein theory. As remarked in this work, 
Wesson [13] and Liko [14] , as well as previous work on this structure[3], have associated 
the source of the 05 field with mass density. A Hamiltonian of the form 

K = ^-g* 1 "^ - eaj(p v - ea u ) (34) 
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can be put into correspondence with K by taking g^ u to be 

r = ^ * (35) 
k + eas 

where k is the common (constant) value of K and K. In this correspondence, the equations 
of notion generated by K through the Hamilton equations, have extra terms, beyond those 
provided by the connection form associated with g^ v ', due to the presence of the gauge fields. 
These additional terms can be identified as belonging to a connection form associated with 
a five dimensional metric, that of a Kaluza-Klein theory. 
We may apply the same procedure to the Hamiltonian 

K = ^9^(P» ~ <MJ(P» ~ <Mv) + ^ (36) 

where g^ v is an Einstein metric, $ is a world scalar field, and W M are to be identified with 
the Bekenstein- Sanders fields for which[6] U U U U = — 1, with = g^ v U v . 

We discuss in Section 4 a class of gauge transformations on the wave functions of the 
underlying quantum theory for which the arise as gauge compensation fields. 

Let us define, as in Eq. (35), the conformally modified metric 



k-$ (37) 

= e" 5 

The "equivalent" Hamiltonian 



K = ^""(Pv ~ <Mj(Pv ~ <Mv) (38) 

then generates, through the Hamilton equations, an equation of motion which corresponds 
to the geodesic equation for an effective Kaluza-Klein metric, as in ref.[8]. 
Now, consider the Hamiltonian 

K K = ^'"'PvPu, (39) 

with the Bekenstein- Sanders metric[7] 

~ gi xv = e -24>( g »» + _ e i<t>um v (40) 

The Hamiltonian Kk then has the form 

K K = e-^g^p^ - 2 sinh 2</>(W> M ) 2 , (41) 

Let us now define a Kaluza-Klein type metric of the form obtained in [7] , arising from 
the equations of motion generated by (38), 

«"-(£ W- (42) 
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Contraction to a bilinear form with the (5D) vectors pa = {px^Ps}, with indices A = v on 
the right and A = \x on the left, one finds 



9 Ab papb = g^P^Pu + 2 P5 (p^) + ( P5 ) 2 g 



(43) 



If we take 




(1 ± y/l - 2# 55 sinh2(/>) , 



(44) 



then the Kaluza-Klein theory coincides with (41), i.e., 



K K = - — 9 PaPb- 



(45) 



As discussed by Wesson [13], Kaluza [15] chose (755 = const, for consistency with electro- 
magnetism, while Wesson [13] makes the more general choice of a world scalar field. In 
particular, the value g 55 = is well defined (as in [8]). 

Since the fields are timelike unit vectors [7], (p^U^) corresponds, in an appropriate 
local frame, to the energy of the particle, close to its mass in the case of a nonrelativistic 
particle, or to the frequency in the case of on-shell photons. It clearly remains to understand 
more deeply the apparently ad hoc choice of p 5 in (44) in terms of a 5D canonical dynamics, 
along with the structure of the 5D Einstein equations for qab that follow from the geometry 
associated with (45). We shall study these questions in a succeeding paper. 

4 The Bekenstein- Sanders Vector Field as a Gauge Field 

Essential features of the Bekenstein- Sanders field [7] of the TeVeS theory are that it 
be a local field, i.e., U^ix), and there is a normalization constraint 



so that the vector is timelike. To preserve the normalization condition (46) under gauge 
transformation, we shall study the construction of a class of gauge transformations which 
essentially moves the U{x) field on a hyperbola with a Lorentz transformation (at the point 



If we think of our underlying quantum structure, which generates the gauge field, as 
a fiber bundle with base x M , then we must think of the transformation acting in such a 
way that the absolute square (norm) of the wave function attached to the base point x M 
preserves its value [9]. 

An analogy can be drawn to the usual Yang-Mills gauge [9] on SU(2), where there 
is a two- valued index for the wave function ip a (x). The gauge transformation in this case 
is a two by two matrix function of x, and acts only on the indices a. The condition of 
invariant absolute square (probability) is 



WU» = -1, 



(46) 



x). 




(47) 
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Generalizing this structure, one can take the indices a to be continuous, so that (47) 
becomes 

J{dU)\ j \dU')U{UM')^{U\x)\ 2 = j \dU)\i){U,x)\ 2 , (48) 

implying that U(U,W) is a unitary operator on a Hilbert space L 2 (dU). Since we are 
assuming that lies on an orbit determined by (48), the measure is 

= (49) 

i.e., a three dimensional Lorentz invariant integration measure. 

Moreover, the Lorentz transformation on is generated by a non-commutative oper- 
ator, and therefore the gauge transformation is non-Abelian. We demonstrate the resulting 
noncommutativity of the operator valued fields, W, after an infinitesimal gauge transfor- 
mation of ths type, explicitly below. 

This construction is somewhat similar to the treatment of the electromagnetic poten- 
tial vector and its time derivative as oscillator variables in the process of second quantiza- 
tion of the radiation field (the energy density of the field is given by these variables in the 
form of an oscillator). 

We now examine the gauge condition: 

( Pfl -eU^U^ = U( Pfl -eU^ (50) 
Identifying Pil with —id/dx^, and cancelling the terms Up^ifj on both sides, we obtain 

W, = UU,V-i - ~U~\ (51) 

in the same form as the Yang-Mills theory [9] . It is evident in the Yang-Mills theory, that 
due to the matrix nature of the second term, the field will be algebra-valued, resulting in 
the usual structure of the Yang-Mills non-Abelian gauge theory. Here, if the transformation 
U is a Lorentz transformation, the numerical valued field would be carried, in the first 
term, to a new value on a hyperbola. However, the second term may well be operator 
valued on L 2 (dU), and thus, as in the Yang-Mills theory, W' M would become nonabelian. 
It follows from (51) that the field strengths 

! ^- 9 ^- S ^ + w -^ ] (52) 

are related to the the field strengths in the transformed form 



according to 

fU x ) = Uf^{x)U~\ (54) 
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just as in the finite dimensional Yang-Mills theories. 
This result follows from writing out, from (51), 

dU' 8U , dU u , OU' 1 

i d 2 U u _ 1 i dU dU- 1 
~ e dx^dx v ~ edx^ dx u ' 

and subtracting the same expression with /i, v reversed. Then add the result to 



QTJ 

r dU , i r dU i dU 

Whenever the combination 
appears, it should be replaced by 

ou- 1 



(56) 



<9a^ 

The result (54) then follows after a little manipulation. 

Now, consider the possibility that this finite gauge transformation leaves U^U^ = — 1. 
We write out 

(lIU^U-i - ^U-^UWU^ - l -^U~i) = -1 - igLu»U-i 
e ox^ e ax^ e ax^ 

i , dU , 

= -i _ « + ^ af/_1 



e cte^ e M dx^ 

^_3U_dU^_ 
e 2 <9a; M <9x M 

(57) 

It may be possible that U can be chosen to make all but the first term in (57) vanish, 
but in the case of finite gauge transformations, it is not so easy to see how to construct 
examples. For the infinitesimal case, it is, however, easy to construct a gauge function 
with the required properties. For 

U = l + iG, (58) 

where G is infinitesimal, (51) becomes 

1 dG 



U^=U,+ i[GM»] + + 0(G 2 ). (59) 
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Then, 



Let us take 



c = -?EKw*).^-^)} 



2 ^ L <9W 7 <9W A / 61 \ 



= -X){wA 7 (W,x),^} 



2 

where symmetrization is required since is a function of U as well as x, and 

This construction is valid in the initially special gauge, which we shall call the "special 
abelian gauge", in which the components of W M commute. The appearance of VC 1 in 
the gauge functions is then admissible since this quantity acts on the wave functions 
< U,x\ip) = ip{U,x) at the point x,in the representation in which the operator W M on 
L 2 (dU) is diagonal. 

Our investigation in the following will be concerned with a study of the infinitesi- 
mal gauge neighborhood of this limit, where the components of W M do not commute, and 
therefore constutite a Yang Mills type field. We shall show in the limit that the corre- 
sponding field equations acquire nonlinear terms, and may therefore suppress the caustic 
singularities found by Contaldi et al [10]. They found that nonlinear terms associated with 
a non-Maxwellian type action, such as (<9 M ZY^) , could avoid this caustic singularity, so 
that the nonlinear terms we find as a residue of the Yang-Mills structure induced by our 
gauge transformation might achieve this effect in a natural way. 

The second term of (60), which is the commutator of G with U^U^ vanishes, since 
this product is Lorentz invariant (the symmetrization in G does not affect this result). 

We now consider the third term in (60). 



= \\N x "< dUXl W + dUXl iV A W (63) 
2 L dx^ dx^ 

dx* dx» J 



There are two terms proportional to 



dx» 

If we take (locally) 



wa 7 (W, x) = u\y(k u x u ), (64) 
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where k v U v = 0, then 

^-U, = k.WJ Xl = 0. (65) 

For the remaining two terms, 



dx^ dx^ 



= N X ^W 



,du Xl du Xl . 



(66) 



+ [U», N^] + zr^LWN* 
1 J dx^ dx» 

Since the commutators contain only terms linear in W M and they have opposite sign, they 
cancel. The remaining terms are zero by the argument (65). The condition U^W 1 = — 1 is 
therefore invariant under this gauge transformation, involving the coefficient u Xl which is a 
function of the projection of rr M onto a hyperplane orthogonal to U^, i.e., a function of 
where kyJA^ = 0. The vector of course, depends on^ (for example, k^ = 14^(14 •b) + b /J ,, 
for some 6 M 7^ 0). 

We now demonstrate explicitly the nonabelian nature of the gauge fields after in- 
finitesinal gauge transformation. With (59), the commutator term in (53) is 

1 3C 1 3C 

MM = (W„ + *[G,W M ] + -q^Mu + i[G,U v ] + 

1 c)C 1 r)C 

- (M „ + i[G , w j + -_)( M , + i[G , w j + -_) (67) 

+ z[W M ,[G,W„]]-z[W,,[G,W M ]], 

where the remaining terms have identically cancelled. Note that this expression does not 
contain any noncommutative quantities. Now, 

[G, U v \ = 2ieu^Uj (68) 

and 



The terms involving [G, U v \ and [G, ZY M ] therefore cancel, so that 

dx v dx^ 



KM = 2%(^-tf) (™) 



We have taken u) X fl = uj x ^(k a x a ), so that 

du xf " 



dx v 
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= (71) 



and therefore 

[Ul,U' v ] = 2i(Ku'\ - k^'\)U x , (72) 

generally not zero. This demonstrates the nonabelian character of the fields. In the Abelian 
limit, we may take u' — > 0, but as we shall a residual nonlinearity, which depends on u" 
may remain in the field equations 

We now consider the derivation of field equations from a Lagrangian constructed with 
the ip's and f^ v f^ v . We take the Lagrangian to be of the form (the indices are raised and 
lowered with g^ v ) 

C = C f + C m , (73) 

where 

£f = -\rf»» (74) 

and 

Q 1 

C m = r(ifr-^(Pv-eU l *)9' tv (Pu-eU v )-*)l> + c.c. (75) 

We shall be working in the infinitesimal neighborhood of the special gauge for Abelian 
Ufx, for which it has the form given in (59) for infinitesinal G. It is therefore not Abelian to 
first order, but we take its variation 8U to be a c-number function, carrying the variation, 
to lowest order, by variation of the first term in (59), and not varying the part of U 
introduced by the infinitesimal gauge transformation (evaluated on the original value of 

u). 

In carrying out the variation of C m , the contributions of varying the ^'s with respect 
to U vanish due to the field equations (Stueckelberg-Schrodinger equation) obtained by 
varying ift* (or ift), and therefore in the variaton with respect to U, only the explicit 
presence of U in (75) need be taken into account. 

Note that for the general case of U generally operator valued, we can write 

r(p» - dA^{p u - tU v )i> = g^((p» - eU^YiPu - d4 v )4>, (76) 

since the Lagrangian density (75) contains an integration over (dU')(dU") (considered in 
lowest order) as well as an integration over (dx) in the action and the operators U are 
Hermitian. In the limit in which U is evaluated in the special Abelian gauge (real valued) , 
and noting that p M is represented by an imaginary differential operator, we can write this 

as 

^VCpa. - MJipv - = -g^ip* + <U»Wip v - eU„)ii, (77) 

i.e., replacing explicitly p^ by —i^d/dx^) = — i<9 M , we have 

SuC m = -^{V>*(<9„ - ieU»)if> - ((^ + ieU^r)^}SU^, (78) 
where we have called g^ v bU v = 5U^, or, 

S u C m = j^U,x)SU^ (79) 
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where j M has the usual form of a gauge invariant current. 

For the calculation of the variation of Cf we note that the commutator term in (52) 
is, in lowest order, a c-number function, as given in (72). 

Calling 

u'\U x EE Vp, (80) 

we compute the variation of 



Then, for 



we compute 



[U' li ,U , v ]=2i(k v v lt -k li v v ) (81) 



SuKM = Su^KMl (82) 

KM = n^v, + K^f) - u)). (83) 



With our choice of k v = U V {U • b) + b v 

dk 



. )u - = 8S(U-b)+U„b\ (84) 



so that 



[W;, U' v \ = 2i{5^{U ■ b) + Uub^v" 
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= O 1 



z.e. 



S u [U , li M = Cr iiV SU r (86) 

The quantity u p is proportional to the derivative of ui*. In the limit that ui,ui' — > 
(cf. (81)), the second derivative, u" which appears in O 1 ^ u may not vanish (somewhat 
analogous to the case in gravitional theory when the connection form vanishes but the 
curvature does not), so that this term can contribute in limit to the special Abelian gauge. 
Returning to the variation of Cf in (74), we see that 

6C f = -d v f^bW + 2if^S[U^,U v }, (87) 

where we have taken into account the fact that \U/j,,U u ] is a commuting function, and 
integrated by parts the derivatives of 5U. With (86) we obtain 

Act 



5C f = -d v j^W + 2ief XlT A \8U» 

Since the coefficient of SW 1 must vanish, we obtain, with (79), the Yang-Mills equations 
for the fields given the source currents 

d u U„ = j»-2ief Xa O x %, (89) 
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which is nonlinear in the fields U^, as we have seen, even in the Abelian limit, where, from 
(78) and (79), 

3„ = - idA^ - + idA^W)^}. (90) 

We point out that this current corresponds to a flow of the matter field; the absolute 
square of the wave functions corresponds to an event density. The coupling e is not 
necessarily charge, and the fields U are not necessarily electromagnetic even in the Abelian 
limit. However, the Hamiltonian (36) leads directly to a Lorentz type force, similar in 
form to that generated by the Hilbert-Einstein action. The dynamics of this system will 
be investigated in a forthcoming paper. 

5. Conclusions 

A map of the type discussed in ref. [8] , of a Hamiltonian containing an Einstein metric, 
generating the connection form of general relativity, and a world scalar field, representing 
a distribution of energy on the spacetime manifold, into a corresponding Hamiltonian with 
a conformal metric (and compatible connection form), can account for the structure of 
the RAQUAL theory of Bekenstein and Milgrom[5]. Furthermore, applying this corre- 
spondence to a Hamiltonian with gauge-type structure, we have shown that one obtains a 
non-compact Kaluza-Klein effective metric which can account for the TeVeS structure of 
Bekenstein, Sanders and Milgrom[7]. 

In order to maintain the constraint condition UpJA^ = — 1 for the Bekenstein-Sanders 
fields, under local gauge transformations, we have introduced a class of gauge of gauge 
transformations on the underlying quantum theory which acts on the Hilbert bundle, 
quite analogous to that arising in the second quantization of the electromagnetic field 
(where the vector potentials and their time derivatives are considered as quantum oscillator 
variables) associated with the values of the gauge fields. The action of this class of gauges 
induces a nonabelian structure on the fields, which therefore satisfy Yang-Mills type field 
equations with source currents associated with matter flow. In the Abelian limit, these 
equations contain residual non-linear terms which may avoid the caustic singularities found 
by Contaldi et al for an electromagnetic type gauge field. 

The phenomenological constraints placed on the TeVeS variables in its astrophysical 
applications and on its MOND limit [16] would, in principle, place constraints on the vector 
and scalar fields appearing in the corresponding Hamiltonian model, for which the additive 
world scalar field corresponds to an energy distribution not associated with electromagnetic 
radiation. 
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